arXiv:1508.01066vl [cond-mat.supr-con] 5 Aug 2015 


Quasiparticle Interference in Fe-based Superconductors Based on a Five-Orbital 

Tight-Binding Model 

Youichi Yamakawa^’0 and Hiroshi Kontani^ 

^Department of Physies, Nagoya University, Furo-eho, Nagoya 464-8602, Japan 

(Dated: August 6, 2015) 

We investigate the quasiparticle interference (QPI) in Fe-based superconductors in both the s++- 
wave and s±-wave superconducting states on the basis of the five-orbital model. In the octet 
model for cuprate superconductors with d^ 2 _y 2 -wsive state, the QPI signal due to the impurity 
scattering ai q = ki — kj (E = |A(fcQ|, z = 1 ~ 8) disappears when the gap functions at ki 
and kj have the same sign. However, we show that this extinction rule does not hold in Fe- 
based superconductors with fully-gapped s-wave state. The reason is that the resonance condition 
E = |A(fci)| is not satisfied under the experimental condition for Fe-based superconductors. We 
perform the detailed numerical study of the QPI signal using the T-matrix approximation, and 
show that the experimentally observed QPI peak around q 2 = (tt, 0) can be explained on the basis 
of both the s++-wave and s±-wave states. Furthermore, we discuss the magnetic field dependence 
of the QPI by considering the Zeeman effect, and find that the field-induced suppression of the peak 
intensity around q 2 can also be explained in terms of both the s++-wave and s±-wave states. 

PACS numbers: 74.70.Xa, 74.20.-z, 74.55.+V 


I. INTRODUCTION 

Since the discovery of Fe-based superconductors^ 
much effort has been devoted to reveal the mechanism 
of high-Tc superconductivity (SC). The mother com¬ 
pounds exhibit structure and antiferromagnetic transi¬ 
tions. These transitions are suppressed by carrier doping 
and then the SC state emerges. In the early theoreti¬ 
cal studies, spin fluctuation mediated s±-wave state, in 
which the SC gap functions change their sign between the 
hole and electron Fermi surfaces (FS), was proposed^^”— 
On the other hand, the orbital fluctuations can induce the 
s++-wave state without sign change in the gap functions 
as discussed in Refs. Figure [1] shows the unfolded 

FS and schematic picture of the (a) 5++-wave and (b) 
5±-wave states. The s++(5±)-wave state is driven by the 
orbital (spin) fluctuations at q 2 = (tt, 0) that corresponds 
to the nesting between hole and electron FSs. 

To distinguish between the s±-wave and s++-wave 
states, various phase sensitive experiments have been per¬ 
formed, such as the impurity effect on Tc^~— the reso¬ 
nant peak by the inelastic neutron scattering^ the co¬ 
herence peak by the nuclear magnetic resonancethe 
quasiparticle interference (QPI) by the scanning tunnel¬ 
ing microscope (STM)^~— and so on. Many theorists 
have preformed theoretical investigations of such experi¬ 
ments based on the realistic five-orbital model. For exam¬ 
ple, the present authors have shown that the robustness 
of Tc against impurities is inconsistent with the 5±-wave 
stateji^ii^ It has been shown that the broad resonant 
peak in the neutron scattering spectrum can be explained 
on the basis of the s++-wave state rather than the 5±- 
wave stateAlso, the absence of the coherence peak at 
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Tc can be explained in terms of both the s++-wave and 
s±-wave states The theoretical study of the QPI signal 
in Fe-based superconductors ware performed by several 
theoretical groups in Refs. [2214^ . 

By using the STM measurement, the information of 
the local density of states can be obtained. The QPI 
signal Z{q^E) is given by the Fourier transformation of 
the tunneling conductance ratio Z{r,E) = 
derived from the STM measurement. The QPI study 
played a crucial role to determine the pairing symme¬ 
try in cuprate superconductorsIn cuprate super¬ 
conductors, the nodal da, 2 _^ 2 -wave SC state is realized. 
There are eight k points {ki : z = I ^ 8) on the FS 
satisfying the relation E = |A(A;i)| for E < A™^^. It 
is called the octet model, and the QPI signal Z{q,E) 
with q = ki — kj emerges due to the impurity scattering 
when A{ki) and A{kj) have the opposite sign, while it 
disappears when A{ki) and A{kj) have the same sign. 
The disappearance of the QPI signal is called the “ex¬ 
tinction rule”. Furthermore, the experimental QPI peak 
is rapidly suppressed by applying a magnetic field. The 
extinction rule and the magnetic field dependence of the 
QPI obtained in cuprate superconductors are well under¬ 
stood in terms of the octet model with dx 2 _y 2 -wdive gap 
symmetry. 

In Fe-based superconductors, many 

experimentali^-— and theoretical^Sri^^ studies 
of the STM have been performed. Hanaguri et al. 
carried out the QPI experiments on Fe(Se,Te) single 
crystal and reported the appearance of a shape peak 
around q 2 = (tt, which is caused by the impurity 

scattering between hole and electron FSs. By analogy 
with the extinction rule in the octet model, the existence 
of the QPI peak around q 2 may indicate that the gap 
functions on the hole and electron FSs have opposite 
sign, i.e., 5±-wave state. Although, many pioneering 
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theoretical studies had been performed for Fe-based 
superconductors, some previous theoretical studies 
assumed over-simplified band structures. Furthermore, 
the QPI signal in the s++-wave state had not been 
studied in detail in previous studies. Therefore, detailed 
theoretical study of the QPI based on a realistic five- 
orbital model in both the s++-wave and s±-wave states 
had been required. 

Hanaguri et al. showed that the intensity of the 
QPI peak around q 2 is slightly suppressed by the mag¬ 
netic field B = 11 T ai E = 1.0 meV^^ However, 
the field-induced change of the QPI peak around q 2 
non-monotonically depends on E] the peak intensity is 
slightly enhanced at E = 0.5 meV and E = 1.9 meV. 
[See Fig. 3S(I) in the Sup plemental Material of Ref. M 
and Fig. 1(A) in Ref. R^.] Therefore, in this paper, we 
discuss the field-induced change of the QPI for wide range 
of E in terms of both the s++-wave and s±-wave states. 

In this paper, we investigate the QPI in Fe-based su¬ 
perconductors on the bases of both the 5++-wave and 
s±-wave states. In the cuprate superconductors with 
dx 2 _y 2 -wave SC state, the QPI signal at q = ki — kj 
disappears when A{ki) and A{kj) have the same sign. 
However, such extinction rule does not hold in Fe-based 
superconductors with fully-gapped s-wave SC state, since 
the resonance condition E = |A(A;i)| = |A(A5j)| is not 
satisfied under the experimental condition E < A™^^ re¬ 
gardless of the sign of the gap functions. We perform the 
detailed numerical study of the QPI signal based on the 
five-orbital model, and find that the experimentally ob¬ 
served QPI peak around q 2 = (tt, 0) appears in both the 
5++-wave and 5±-wave states. Furthermore, we discuss 
the magnetic field dependence of the QPI by considering 
the Zeeman effect, and find that the field-induced change 
of the peak intensity around q 2 can also be explained 
in terms of both the s++-wave and 5±-wave states. In 
conclusion, it is difficult to distinguish between the s++- 
wave and s±-wave states from the QPI experiments in 
Fe-based superconductors. 


II. FORMULATION 


A. Quasiparticle Interference 


The tunneling conductance dI/dV{r^ V) at position r 
and voltage V is approximately proportional to the lo¬ 
cal density of states p{r,E) at energy E = namely, 
dlldV{r^V) (X \M{r)f p{r,E), where we set the unit of 
charge e as one. M(r) is the tunneling matrix element 
between the sample surface and the STM tip. In the pres¬ 
ence of the impurities, we can drop the factor M(r) and 
obtain the information of the density of states by taking 
the ratio Z{r,E) between the conductance measured at 
-\-V and — U as follows 


Z{r,E) 


dI/dV{r,EV) _ p{r,EE) 
dI/dV{r,-V) ~ p{r, -E) 
p%+E) r 5p{r,+E) 6p{r,-E) 

p%-E)[ ^ pO{r,+E) p^ir,-E) 


,( 1 ) 


where p^{E) is the averaged density of states and 6p{r, E) 
describes the spatial modulation defined as 6p{r,E) = 
p{r,E) — p^{E). The Fourier transformed conductance 
ratio is called the QPI signal, which is given by 


Z{q,E) 


P°{+E) 

p^{-E) 


V(g,+-E) 

p%+E) 


6p{q,-Ey 


( 2 ) 


where q is a scattering wave vector. When the system 
is uniform, Z{q^E) is zero except for gr = 0. We can 
obtain the information on the SC gap symmetry since 
the momentum dependence of Z{q^E) reflects the sign 
of the gap functions. 


B. Model Hamiltonian and Green Function 



FIG. 1. (Color online) Fermi surfaces and gap structures in 
the (a) s++-wave and (b) s±-wave states. The arrows denote 
scattering wave vectors. The scattering vector q 2 ~ (7r,0), 
which is equal to the nesting vector, connects hole and elec¬ 
tron Fermi pockets. qi ~ (0,0) corresponds to the intra¬ 
band scattering, and qs ~ (tt, tt) corresponds to the scattering 
within electron or hole Fermi pockets. 


The five-orbital tight-binding Hamiltonian is given by 

E (3) 

where ^ ^ {ck,i,a) is the creation (annihilation) operator 
of a Fe 3d electron with wave vector k, orbital I and 
spin a. is given by the Fourier transformation of 

the hopping integrals introduced in Ref. Q. The energy 
dispersion of band b is obtained as a eigenvalue of 
by unitary transformation, 

1 , 1 ' 

where Uh,i,b is an element of the unitary matrix obtained 
as the eigenvector. The obtained Fermi Surface is shown 
in Fig. [TJ 
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Now we study the SC state. In a single-orbital model, 
the BCS Hamiltonian is simply given by 

H = (5) 

k 

where 

Here, we define the Pauli matrices and which act 
in particle-hole space and spin space, respectively. For 
example. 
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(7) 


Then, the Nambu Hamiltonian Hk for a single-orbital 
model is given by 


Hh = efcfa - Hds + A^fi, 


(8) 


where B is the Zeeman splitting energy by the magnetic 
field and is the singlet gap function. 

In the five-orbital model, the Nambu Hamiltonian is 
written as follows: 

Hk = Hlh - , (9) 


where is the 5x5 unit matrix in the orbital space. In 
the case of the five-orbital model, the Nambu Hamilto¬ 
nian is given by the 20 x 20 matrix form. A^^^ is 5 x 5 
matrix form singlet gap function in the orbital space, and 
its matrix element is obtained by the unitary transforma¬ 
tion of the band-basis gap function Ak,b as 

( 10 ) 

b 

Then, the Green function in the clean limit is given by 

Glico) = (uj - HkY" 

= {lo- H°fs + BE5&3 - , ( 11 ) 

and the local density of states without randomness is 
given by 




T 0 +T 3 


Glico) 




( 12 ) 


where 7 is the quasiparticle damping rate. 

When we consider the impurity scattering, the Green 
function is obtained by using the T-matrix approxima¬ 
tion as follows: 


Gk,k'{^) — G^{^)^k,k' + ^Gk,k'{^)^ (13) 


where 


5Gk,k'H = Gl{Lo)fk,k'{Lo)Gl,{Lo). (14) 

For a single impurity, the T-matrix is obtained by solving 
the following self-consistent equation, 

Tk,k'{^) = Ik,k' + ^ Ik,k” G%^ ,k'{^), (15) 

k” 


where Ik,k' is the impurity potential of a single impurity. 
The modulation of the density of states induced by the 
impurity scattering is given by^^^ 


dp{q,uj) = - E (5Gfe,fe+q(w) 

k 


(16) 


where G is represented in the orbital basis. This treat¬ 
ment is exact for the case of low impurity concentration 
< 1 . 

In this paper, we consider the non-magnetic impurity 
since the QPI due to the magnetic impurity scattering is 
subdominant for T = 0^^ According to the band calcula¬ 
tions, the impurity potential in Fe-based superconductors 
is screened and well-localized^^ That is, the impurity 
scattering matrix in the orbital space is /^-independent. 
When the Fe-site substitution is considered, the impurity 
potential is given as 

/imP = /imp473^ (17) 

and then, the T-matrix becomes /^-independent and it is 
simply given by 

f (w) = (l - /™P5°(w)) J™P, (18) 

where g^{uj) = ^ G^{uj) is the local Green function 
in the 20 x 20 matrix form. 


III. RESULT 


A. Simple Analytical Calculation 


In this section, we analytically show that the extinction 
rule, which tells that the non-magnetic impurity scatter¬ 
ing between FSs with same sign gap functions does not 
contribute to the QPI, does not hold in fully gapped 5 - 
wave SG state. 

Here, we verify the case with the particle-hole sym¬ 
metry p^(+T') = p^{—E). Then, Z{q^E) in Eq. (|2]) is 
simplified as 


Z{q^0,E) = 


26p^^^{q,E) 

P°{E) ’ 


( 19 ) 


where 


E) = +E) Sp{q, E) ^ 
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When we consider the scattering due to non-magnetic 
impurities with a weak scalar potential the T- 

matrix is given by T From Eq. ([16]) . the 

modulation of the density of states for T ^ is 

given by 


5p-<^\q,E) 




ttAT 


^impjimp 

2ttN 


k,b,b' 


E‘^ + 


^k,b^k-\-q,b' 


— Afc /.A 


fc, 6 ^fc+q ,6 


As shown in Fig. [U the other QPI signal can arise 
around qi = ( 0 , 0 ) due to intra-band scattering, and 
around gs = (tt, tt) due to inter-band scattering between 
hole-FSs or electron-FSs. Experimentally, the QPI peak 
around gs is enhanced by the external magnetic field. 
However, both the QPI peaks around qi and are 
caused by the scattering between hole-pockets and be- 
E=E+i^ tween electron-pockets. These QPI peaks are not useful 
for the purpose of distinguishing between the s++-wave 
and s±-wave states. 




E 

I 




E=E-\-i^ 


( 21 ) 


B. QPI for the Weak Impurity Potential Case 


In the last line, we utilized the functional form of the 
Green function in the band-diagonal basis, and ^ = 
eL + A 2 ^ is the energy of a quasiparticle in band b. In 
Eq. ([ 2 T]) . the main contribution originates from the case 


that both k and k q are on FSs = ^k-\-q,b' 
In this case, the contribution is simplified as 


6p°<^^{q,E)^- 


E 

k ,k-\-q^FS 


Im 


AfcAfc-|-q 


(E^ - Al){E^ - 

In cuprate superconductors, the nodal d. 


= 0). 


( 22 ) 

E=E-\-ij 


-y^ 


SC state is realized. Under the experimental condition 
E < only the eight k points : i = 1 ^ 8 ) satisfy 

the relation E = |AfcJ. It is called the octet model^“— 
and ki ks are shown in Fig. n a). Sp^'^^{q,E) can 
be very large for q = ki — kj since the denominator in 
Eq. (| 2 ^ is almost zero for k ^ ki. On the other hand, 
the numerator is sensitive to the sign of the gap func¬ 
tions: the numerator has finite value 2E^ when the gap 
functions A^ and A^+g have opposite sign, but it be¬ 
comes zero for the same sign case. Therefore, the QPI 
peak disappears when the gap functions at k and k q 
have the same sign, which is called the extinction rule. 

In contrast, such extinction rule does not hold in the 
fully gapped s-wave SC state realized in Fe-based super¬ 
conductors, under the experimental condition E < A™^^. 
We focus on the QPI peak around q2 = (tt, 0) which cor¬ 
responds to the inter-band scattering between the hole 
and electron FSs. Using the gap functions on hole FS 
Ah and electron FS Ae, is given by 


5p^^^{q2,E) (X -Im 


E‘^ - AhAe 


{E^ - Al){E^ - Al) 


(23) 


E=E-\-ij 


In the case of F < |Ah,e |5 both the numerator and de¬ 
nominator have finite value regardless of the signs of Ah 
and Ae- That is, Sp'^^^{q 2 ^ E) is finite even for AhAe > 0. 
Therefore, the extinction rule does not hold in Fe-based 
superconductors, and the QPI peak around q 2 is ex¬ 
pected to appear in both the s++-wave and s±-wave 
states. We will numerically verify the violation of the 
extinction rule for the q 2 signal by analyzing the five- 
orbital model in later sections. 


In this and subsequent sections, we numerically calcu¬ 
late the QPI signal using Eq. m- Here, we discuss the 
QPI due to a weak impurity potential = 0.1 eV, and 
show that the QPI peak around q 2 is actually obtained 
in both the 5 ++-wave and s±-wave states for various pa¬ 
rameters. Hereafter, we set Aq = 0.02 eV, = 0.01, 
7 = Ao/4 and N = 256 x 256. We confirmed that the ob¬ 
tained results do not change qualitatively for 7 = Ao/8. 

Figure [2] shows the intensity map of the QPI, |E(gr, F)|, 
at zero field. First, we discuss the (i) isotropic single-gap 
case with |Ah| = |Ae| = Aq: Figures [2fa) and (b) show 
the results obtained in the 5++-wave and 5±-wave states, 
respectively. Considering the experimental condition in 
Ref. UH, we set E = Ao/2. In the 5++-wave state (a), the 
sharp QPI peak around q 2 clearly appears as expected 
from Eq. ([23]) . Therefore, the extinction rule does not 
hold in Fe-based superconductors. In the s±-wave state 
(b), the strong QPI peak accompanied by the large halo 
structure is obtained around q 2 . That is, it is difficult to 
distinguish between the s++-wave and 5±-wave states by 
the presence or absence of the QPI peak around q 2 . 

In reality, | Ah | and | Ae | are different in usual Fe-based 
superconductors. For example, ^ 2 is re¬ 

ported in electron- and hole-doped BaFe 2 As 2 i^^i^^ In the 
Fe(Se,Te) sample used in the QPI experiments^ the re¬ 
lations Amin ^ 1 nieV and Amax = 2^4 are expected 
from the tunneling conductance measurement. There¬ 
fore, we show the results for the (ii) isotropic two-gap 
case with |Ah| = 2Ao and |Ae| = Aq in Figs. [2](c) and 
(d). In this case, there is no large difference from the 
single-gap case shown in Figs. n a) and (b). Similar re¬ 
sults are obtained when |Ah| = Aq and |Ae| = 2Ao. 

In Figs. m e) and (f), we also show the (iii) strongly 
anisotropic gap case with |Ah| = Aq and |Ae| = (1 ± 
cos2^)Ao. Anisotropic-gap functions are reported on a 
hole-FS in heavily K-doped BaFe 2 As 2 ^ and on the elec¬ 
tron FSs in some Fe(Se,Te) systems j^^d^ In this case, the 
peak around q 2 exists and its shape in the s++-wave state 
becomes similar to the one in the s±-wave state. There¬ 
fore, it is difficult to distinguish between the s++-wave 
and s±-wave states from the existence of the QPI signal 
around q 2 . 

Experimentally, the QPI peak intensity around q 2 is 
slightly suppressed by the magnetic field B for E = 
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^++-wave state ^+-wave state 





FIG. 2. (Color online) Intensity map of the QPI \Z{q,E)\ 
under zero field S = 0atF^ = Ao/2 due to the non-magnetic 
impurity scattering with potential = 0.1 eV. The left 
and right panels show the results in the s++-wave and s±- 
wave states, respectively. (a),(b) Case (i): Isotropic single¬ 
gap case with Ah = ±Ae = Aq. (c),(d) Case (ii): Isotropic 
two-gap case with Ah = 2Ao and Ae = ±Ao. (e),(f) Case 
(iii): Strongly anisotropic gap case with Ah = Aq, Ae = 
±(1 + cos 2^) Ao around k = (tt, 0) and Ae = ±(1 — cos 2^) Aq 
around k = (0,7r). The vertical broken lines represent the 
path of the linecuts in Fig. [3l 


1.0 meV < Here, we discuss B— and E— de¬ 

pendencies of \Z{q,E)\ in detail, and show that the ex¬ 
perimental suppression of the q 2 peak for E ^ A™^^ can 
be explained in both the s++-wave and s±-wave states. 
Previously, two kinds of the field-induced suppression ef¬ 
fects have been discussed by Coleman et (A) Im¬ 

purities are masked by vortices under the magnetic field, 
and then the impurity scattering rate is reduced. Also, 
(B) the Zeeman effect changes the electronic state and 
modifies the impurity scattering. The former mechanism 
would suppress the QPI intensity around q 2 regardless of 
the sign of the gap functions. However, in the QPI ex¬ 
periments for Fe(Se,Te) in Ref. M, it was reported that 
the effect (B) would be dominant, since the field-induced 
changes are almost spatially uniform. Therefore, in this 
paper, we study only the effect (B). 

Figures [3l^a) and (b) show the \Z{q^E)\ in the single¬ 
gap case with |Ah| = |Ae| = Aq [case (i)] from q = 
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FIG. 3. (Color online) Linecuts from \Z{q, E) \ map for = 
0.1 eV. The solid and dotted lines represent \Z(q, E)\b=o and 
\Z{q,E)\B =Ao/ 2 , respectively. The left and right panels show 
the results in the s++-wave and s±-wave states, respectively. 
The path is shown in Fig.O (a),(b) Case (i): Isotropic single¬ 
gap case with Ah = ±Ae = Aq. (c),(d) Case (ii): Isotropic 
two-gap case with Ah = 2Ao and Ae = ±Ao. (e),(f) Case 
(iii): Strongly anisotropic gap case with Ah = Aq, Ae = 
±(1 + cos 2^) Ao around k = (tt, 0) and Ae = ±(1 — cos 2^) Aq 
around k = (0,7r). The curves in all the hgures are vertically 
shifted to make them visible. 


(tt, —tt) to (7r,7r). The path is shown in Fig. [2] by the 
vertical dashed lines. The solid and dotted lines repre¬ 
sent the results for H = 0 and B = Ao/2, respectively. 
In the (a) s++-wave state, the QPI peak around q 2 is 
not sensitive to B and E. On the other hand, in the 
(b) s±-wave state, the q 2 peak is drastically suppressed 
by B. Figures m c) and (d) show the results obtained 
for the two-gap case [case (ii)]. In this case, the QPI 
peak around q 2 is suppressed by H for E" ^ Aq in both 
the 5++-wave and 5±-wave states. However, the field 
suppression of the QPI peaks is much larger in the 5++- 
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wave state. Figures [3l^e) and (f) show the results for the 
strongly anisotropic gap case [case (iii)]. In this case, 
\Z{q,E)\ in the s++-wave state shows very complex B 
dependence. 

In summary, in the s±-wave state, the QPI peak 
around q 2 is clearly suppressed in all cases (i)-(iii). 
In the s++-wave state, this peak intensity is also sup¬ 
pressed in the two-gap case (ii). Therefore, the field- 
induced suppression of the QPI around q 2 can be ex¬ 
plained in terms of both the 5 ++-wave and s±-wave 
states. Experimentally, the SC gaps are fully opened in 
the Fe(Se,Te) sample used for the QPI experiments, and 
relation ^ ^mm ^ ^ meV is expected, since the 

estimated value of 2 A™^^/Tc < 2 is much smaller than 
the BCS value 3.53. In addition, the tunneling conduc¬ 
tance has the sharp gap edge peak at E ^ 1.7 mV and an 
additional peak at about 4 mV. If the latter peak arises 
from the SC gap, A^^^/A™^^ > 2 is expected. There¬ 
fore, the isotropic two-gap case with |Ah| = 2|Ae| [case 
(ii)] would correspond to Fe(Se,Te). 


C. QPI for the Strong Impurity Potential Case 

In this section, we consider the QPI due to a strong im¬ 
purity potential = 1 eV, which corresponds to Fe- 

site substitution. Since the residual resistivity takes the 
maximum for ^ +1 eV, P^^ = +1 eV corresponds 
to the unitary limit in Fe-based superconductors 
Here, we show the result only for the isotropic two-gap 
case with |Ah| = 2 Ao and |Ae| = Aq [case (ii)]. 

FigureslUa) and (b) show the \Z{q, E) \ map for = 
— 1 eV in the case of the s++-wave and s±-wave states, 
respectively. Also, Figs. |4l[c) and (d) show the ones for 
Jimp ^ j; ^ and B = 0. The 

obtained QPI map is qualitatively similar to the ones in 
the weak potential case shown in Fig. [ 2 l and the QPI 
peak around q 2 appears in both the s++-wave and s±- 
wave states. Therefore, the extinction rule does not hold 
in Fe-based superconductors regardless of the magnitude 
of the impurity potential. 

Figures [SKa) and (b) show \Z{q,E)\ from q = (tt, —tt) 
to ( 7 r, 7 r) for P^p = — 1 eV. The solid and dotted lines 
represent the results for P = 0 and Ao/2, respectively. 
For P ^ Aq, the QPI peak around q 2 is suppressed by 
the magnetic field in both the s++-wave and 5 ±-wave 
states. 

FiguresISJc) and (d) show the results for P^p = -hi eV. 
In the s++-wave state, the QPI peak around q 2 is insensi¬ 
tive to B and E. On the other hand, in the <s±-wave state, 
the QPI signal shows very strong E dependence, and the 
QPI intensity becomes very small for E > O.SAq even for 
P = 0. However, such behaviors have not been observed 
experimentally. As results, in both the s++-wave and 
s±-wdive states, the obtained results for P^^ = +1 eV 
are inconsistent with experimentsTherefore, impu¬ 
rities with weak potential will be responsible for the QPI 
signal in Fe(Se,Te). 
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FIG. 4. (Color online) Intensity map of\Z{q, E) \ at P = Ao/2 
and P = 0 in the isotropic two-gap case with Ah = 2Ao 
and Ae = ±Ao [case (ii)]. The left and right panels show 
the results in the s++-wave and s±-wave states, respectively. 
(a),(b) P""P = -1 eV. (c),(d) P™P = +1 eV. The vertical 
broken lines represent the path of the linecuts in Fig. [5l 


In the above discussion, we have ignored the change 
of Pc due to the impurity scattering. We have shown 
that the s±-wave state with the original SC transition 
temperature Pco = 30 K is completely suppressed when 
the residual resistivity reaches ^ where 

z~^ = m*/m is the mass-enhancement factor due to the 
self-energy. When = +1 eV, the residual resistiv¬ 
ity for = 0.01 is about 20/if2cm in Fe-based super¬ 
conductors. Therefore, the 5 ±-wave state is very fragile 
against impurity. 


IV. DISCUSSION 


A. Violation of the Extinction Rule 


As shown in Sec. imi the QPI peak around q 2 is re¬ 
alized even in the 5 ++-wave state. The reason is that 
the numerator in Eq. ([23]) has finite value under the 
experimental condition P < |Ah,e|- Thus, the extinc¬ 
tion rule in the octet model for cuprate superconductors 
(|A(A;i)| = P < = 1^8), which tells that the 

QPI signal at gr = — kj disappears if A{ki) = A(/cj), 

does not hold in Fe-based superconductors under the ex¬ 
perimental condition P < Aq. As shown in Fig. jS] the 
QPI signal around q 2 still exists even at P = Aq in the 
s++-wave state due to the finite quasiparticle damping 
7 . For these reasons, we can not distinguish between 
the s++-wave and s±-wave states from the presence or 
absence of the QPI peak around q 2 . 
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FIG. 5. (Color online) Linecut from \Z{q^ E)\b=q (solid lines) 
and \Z{q^ E)\b=Aq /2 (dotted lines) maps in the isotropic two- 
gap case with Ah = 2Ao and Ae = =bAo [case (ii)]. The path 
is shown in Fig.[4l The left and right panels show the results in 
the s++-wave and s±-wave states, respectively. (a),(b) = 

— 1 eV. (c),(d) = +1 eV. The curves in all the figures are 

vertically shifted to make them visible. 


the s++-wave state maybe due to the very large difference 
in the band structure. 


V. SUMMARY 

In summary, we investigated the QPI in Fe-based su¬ 
perconductors in both the s++-wave and s±-wave states. 
In the octet model (|A(A;i)| = E < = 1^8) for 

cuprate superconductors with (ia, 2 _^ 2 -wave SC state, the 
QPI signal around q = ki — kj disappears when A{ki) 
and A{kj) have the same sign. However, this extinction 
rule is not hold in Fe-based superconductors with fully- 
gapped s-wave SC state. The reason is that the resonance 
condition, in which the denominator of the integrand in 
Eq. (|22j) becomes zero at some /c, does not satisfied under 
the experimental condition E < |Ae,h|- We performed 
the detailed numerical study of the QPI signal on the 
basis of the five-orbital model and found that the exper¬ 
imentally observed QPI peak around q 2 = (tt, 0) can be 
explained in terms of both the s++-wave and 5±-wave 
states. Furthermore, we discussed the magnetic field 
dependence of the QPI by considering the Zeeman ef¬ 
fect, and found that the suppression of the peak intensity 
around q 2 by the magnetic field can also be explained in 
terms of both the s++-wave and 5±-wave states. There¬ 
fore, it is difficult to distinguish between the 5++-wave 
and s±-wave states from the QPI experimental date for 
Fe-based superconductors. 


B. Comparison with Previous Studies 

In Ref. [ 23 , Sykora and Coleman investigated the QPI 
in the 5±-wave state by using a two-band model. They 
showed that the QPI peak around q 2 emerges for H = 0 
due to the non-magnetic impurity scattering in the weak 
potential limit, and its intensity is suppressed by the Zee- 
man effect under the magnetic field B = Aq. It is con¬ 
sistent with the result of the present study for the weak 
potential case based on the five-orbital model. Also, to 
analyze the unitary scattering case, they phenomenolog¬ 
ically treated the resonant scattering due to the multi¬ 
ple scattering process, and proposed that the QPI signal 
around q ^, = (tt, tt) is enhanced by B due to the reso¬ 
nant scattering. However, we cannot obtain such behav¬ 
ior in the present study using T-matrix approximation 
for = +1 eV. 

In Ref. [ 2 ^, Gao et al discussed the magnetic field de¬ 
pendence of the QPI due to the vortex, which is not con¬ 
sidered in the present study. Interestingly, they showed 
that the strong and sharp QPI peak around qs is caused 
in both the s++-wave and s-t-wave states by the Andreev 
scattering due to the vortices. Experimentally, however, 
the field-induced change is almost spatially uniform, indi¬ 
cating that the impurity scattering is more important 
In Ref. [ 2 ^, the QPI peak around q 2 was not obtained in 
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Appendix A: QPI in Cuprate Superconductors 

In the QPI measurement for the cuprate by Hanaguri 
et al^ it was shown that the QPI signals due to the im¬ 
purity scattering between k points with opposite sign gap 
functions are strongly suppressed by the magnetic field. 
Since the suppression in the “matrix region” (far from 
vortex) is stronger than the one in the “vortex region” 
(near the vortex core), the Zeeman effect would be impor¬ 
tant. In this appendix, we investigate the magnetic field 
dependence of the QPI in cuprate superconductors with 
nodal dx2_y2-w8ive SC state, A^ = Ao{coskx — cosky)/2, 
using the T-matrix approximation in the case of weak 
impurity potential = 0.1 eV, and show that the ex¬ 
perimentally observed suppression can be explained by 
the Zeeman splitting scenario. 

Eigure [6][a) shows the ES and the gap function in 
cuprate superconductors. The eight wave vectors ki {i = 
1 ^ 8) on the ES satisfy the relation E = |AfcJ < Aq. 
The scattering vectors ^^1,4,5 (^2,3,6,?) connect the two 
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FIG. 6. (Color online) (a) FS of cuprate and the d^ 2 _y 2 - 
wave SC gap. The wave vector (i = 1 ~ 8) satisfies the 
relation E = | A^. |, and qi is scattering vector. < 3 ^ 2 , 3 ,6 ,7 (<?i, 4 , 5 ) 
connect the fc-points on FS with opposite (same) sign gap 
functions, (b) Intensity map of \Z(q,E)\ ai E = Ao/2 and 
B = 0. The scattering vector qi is equivalent to qi. (c) 
\Z{q,E)\ for B = Ao/8. (d) Field-induced change given by 
\^Iq, E)\b=Ao/ 8 — \Z{q, E)\b=o- 


FIG. 7. (Color online) (a)(b) Intensity map of the QPI 
\Z(q, E)\b=o due to the band represented weak impurity po¬ 
tential = 0.1 eV in the (a) s++-wave and (b) s±-wave 

states. (c)(d) The magnetic field-induced change of the QPI 
signal given by \Z{q, E)\b=Ao/ 2 - \Z{q,E)\o in the (c) s++- 
wave and (d) s±-wave states. We set E = 0.7Ao, Ah = 2Ao 
and Ae = ±Ao. 


/c-points with same (opposite) sign gap functions. Ex¬ 
perimentally, the QPI signals are obtained at ^^ 2 , 3 , 6 ,? foi* 
zero field, and they are suppressed by applying a mag¬ 
netic field^^ Figure [U^b) shows the numerical results of 
the QPI intensity map \Z{q^ E)\b=o without magnetic 
field. We use the parameters given in Ref. [3l|. The strong 
QPI peaks appear at q 2 , 3 , 6 , 7 - Figures [6fc) and (d) show 
the QPI with magnetic field, \Z{q, E)\b=Ao/ 8 ^ field- 
induced change given by \Z{q, E^)|b=Ao/8 “ ^)|b=o, 

respectively. In this case, the QPI signal shows re¬ 
markable field dependence and its peaks at q 2 , 3 , 6,7 are 
strongly suppressed by the Zeeman effect. This result 
is consistent with the experimental results for cuprate 
super conductor s 

Appendix B: QPI due to Simplified Impurity 
Potential 

In the above discussion, we have investigated the QPI 
due to the orbital diagonal impurity potential in Eq. (13. 
In this case, the impurity potential has complex k- 
dependence in the band basis. In this appendix, we con¬ 
sider the QPI due to a simple constant impurity potential 
in the band basis, 

/■band _ j ^ ^ ) /-p-i \ 

~\r {b^ b') ’ 

where b = b' and b ^ b' terms correspond to intra¬ 
band and interband scattering, respectively. Hereafter, 


we study the QPI in the weak potential case with I = 
r = 0.1 eV. 

Eigures [Tl^a) and (b) show the QPI intensity map in 
the s++-wave and s±-wave states, respectively. We set 
E = 0.7Ao, |Ah| = 2Ao, and |Ae| = Aq. The QPI 
peak around q 2 appears in both the 5++-wave and s±- 
wave states. Eigures [TKc) and (d) show the field-induced 
change \Z{q, E)\b=Ao /2 - \Z{q,E)\B=o in the s++-wave 
and s±-wave states, respectively. The obtained results 
are qualitatively consistent with the orbital diagonal po¬ 
tential case shown in Eigs. [2fc), [2fd), [3jc), and [3fd) in 
the main text. Therefore, in the weak potential case, 
the obtained QPI signal is insensitive to the nature of 
impurity potential. 

However, the impurity potential in Eq. m gives an 
erroneous result in the unitary regime / ^ 00 , that is, 
the T-matrix becomes band diagonal except for 

I'/I = 1. Due to this model artifact, the QPI peak 
around q 2 disappears in the unitary limit. Eor the same 
reason, Tc in the s±-wave state is almost unchanged by 
impurities in the unitary regime lp^{0) > How¬ 

ever, such erroneous model artifact is revised by using a 
realistic potential in Eq. That is, the QPI peak 

around q 2 appears and Tc in the s±-wave state is fragile 
against impurity even in the unitary regime. 
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Appendix C: Another Two-Gap Case with 

|Ah| = 1.5|Ae| 

In the main text, we discussed the field-induced sup¬ 
pression of the QPI peak intensity around q2 in the 
isotropic two-gap case with |Ah| = 2|Ae|. Here, we 
show another two-gap case with |Ah| = 1.5|Ae|. The 
obtained results are qualitatively the same as the results 
for I Ah I = 2|Ae| in the main text. 

Figure [8] shows the \Z{q^E) \ from gr = ( tt , — tt ) to (7r,7r) 
for H = 0 (solid lines) and B = Ao/2 (dotted lines). 
In the (a) s++-wave and (b) s±-w8ive states, the QPI 
intensity for = 0.1 eV around q2 is suppressed by B 
for E Aq. 

Figure [Hfc) shows the \Z{q,E)\ in the s++-wave state 
for P^^ = —1 eV. In this case, the QPI intensity at 
just q 2 is strongly enhanced by H at ^ Aq, whereas 
the integrated intensity around q2 is suppressed. Such 
field-induced enhancement at just q 2 for £" = Aq is not 
universal since the q 2 peak is suppressed by H for | Ah | = 
2|Ae| as shown in Fig. [5fa) in the main text. However, 
the obtained field-induced enhancement at just q2 may be 
consistent with the experimental result. Experimentally, 
the QPI signal for E = 1.0 meV is suppressed by B 
around q2^ but a slight enhancement is observed at just 
q2 as shown in Fig. 1(A) in Ref. [l^ . 

Figure [8] (d) shows the \Z{q,E)\ in the s±-wave state 
for = —1 eV. Also, Figs, [^e) and (f) show the 

ones for P^p = +1 eV. In all cases (d)-(f) in Fig. [51 the 
obtained results are almost same as the cases (b)-(d) in 
Fig. [5] in the main text for |Ah| = 2|Ae|. 

Therefore, the obtained results for |Ah| = 1.5|Ae| are 
qualitatively same as the ones for |Ah| = 2|Ae| in the 
main text. The field-induced enhancement at just q2 
for /^™P = — 1 eV in Fig. |8l^c) may be consistent with 
experimental result, although it is sensitive to model pa¬ 
rameters. 



FIG. 8. (Color online) \Z(q,E)\ along q — {7T,qy) in the 
isotropic two-gap case with Ah = 1.5Ao and Ae = ±Ao. 
The solid and dotted lines represent P = 0 and B = Ao/2, 
respectively. The left and right panels show the results in the 
s++-wave and s±-wave states, respectively. (a),(b) P^^ = 0.1 
eV. (c),(d) P^^ = -1 eV. (e),(f) P^^ = -hi eV. The curves 
in all the figures are vertically shifted to make them visible. 
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